Interrelation between the pseudo-gap and the incoherent quasi-particle features of 

high-T c superconductors 
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Using a scenario of a hybridized mixture of localized bipolarons and conduction electrons, we 
demonstrate for the latter the simultaneous appearance of a pseudogap and of strong incoherent 
contributions to their quasi-particle spectrum which arise from phonon shake-off effects. This can 
be traced back to temporarily fluctuating local lattice deformations, giving rise to a double-peak 
structure in the pair distribution function, which should be a key feature in testing the origin of these 
incoherent contributions, recently seen in angle resolved photoemission spectroscopy (ARPES). 
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The appearance of a pseudogap, accompanied by a 
predominantly incoherent quasi-particle spectrum in cer- 
tain parts of the Brillouin zone fjj], is considered to be 
amongst the most significant signatures of high-T c super- 
conductors (HT C SC) which may contain the key of our 
understanding of these materials. As suggested earlier, 
the large incoherent part of the quasi-particle spectrum 
might come from a coupling of the electrons to collec- 
tive modes such as spin fluctuations 0. We shall dis- 
cuss and defend in this Letter a similar point of view 
based on a scenario of a mixture of intrinsically localized 
bipolarons and coexisting itinerant electrons, hybridized 
with each other via charge exchange, permitting bipo- 
larons to disintegrate into pairs of conduction electrons 
and vice-versa to reconstitute themselves in an inverse 
process. The location of the bipolarons in high-T c mate- 
rials might be sought in the highly polarizable dielectric 
layers adjacent to the CuC>2 planes or possibly inside the 
polaronic stripes || in those planes themselves - the re- 
mainder of those CuOi planes forming the subsystem 
housing the itinerant electrons. Taking the bipolarons as 
quasi-particles without any internal structure, such a sce- 
nario is described by the so called Boson-Fermion model 
(BFM) which has led us to the prediction of pseudo- 
gap features in the quasi-particle spectrum driven 
by strong local electron-pair correlations. In the present 
Letter we extend our previous studies by taking into ac- 
count the internal polaronic structure of the bipolaronic 
Bosons, as being composed of charge and lattice vibra- 
tional degrees of freedom, locked together in a coherent 
quantum state. A bipolaronic Boson localized on a site i 
is represented by 

b+e-« a *-<)\0)\Z(X)) = b+\0)\*(X-Xo)). (1) 

where the phonon operators a\ +%> correspond to local lat- 
tice deformations. The hard-core Bose operators b\ + ' 
describe pairs of electrons which are self-trapped inside 



locally deformed clusters of atoms, characterized by de- 
formed harmonic oscillator states \&(X — Xq)) with equi- 
librium positions shifted by Xq = 2ay/h/2Mu> (lu de- 
notes the characteristic frequency and M the mass of the 
oscillators). The strength of the coupling of the charge 
carriers to local lattice deformations, ultimately leading 
to bipolaron formation, is given by fioj^a. Such physics is 
described in terms of the following generalization of the 
original BFM: 

+ (As - 2/x) Kbi + + c+c+h] 



(2) 



Here are Fermionic operators referring to itiner- 

ant electrons with spin a. The bare hopping integral 
for the electrons is given by t, the bare Fermionic half 
band width by D, the Boson energy level by As and 
the Boson-Fermion pair-exchange coupling constant by 
v. The chemical potential fi is common to Fermions 
and Bosons. The indices i denote effective sites involving 
molecular units made out of adjacent molecular clusters 
of the metallic Fermionic and dielectric Bosonic subsys- 
tems. Because of the small overlap of the oscillator wave 
functions at different sites we may, to within a first ap- 
proximation, consider the Boson and Fermion operators 
as commuting with each other. 

The original BFM, given by the first two lines in 
Eq.(2), has been investigated in great detail as far as 
the opening of the pseudogap is concerned and as far 
as this affects the thermodynamic, transport and mag- 
netic properties H . The opening of the pseudogap in the 
Fermionic density of states was shown to be driven by 
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the onset of local electron pairing without any supercon- 
ducting long range order. Even without treating the gen- 
eralized BFM within the self-consistent conserving ap- 
proximation, used in those studies of the original BFM , 
we find that the atomic limit of this generalized BFM 
already gives us clear indications on the interrelation be- 
tween the opening of the pseudogap and the appearance 
of predominantly incoherent quasi-particle features, as 
seen in ARPES studies. 

In order to set the scale of the various parameters in 
this model we measure them in units of D, which for typ- 
ical HT C SC is of the order of 0.5 eV. As in our previous 
calculations of the original BFM, we choose v such that 
the pseudogap opens up at temperatures of the order of 
a hundred degrees K. We take v = 0.25 for the present 
study. We furthermore choose a = 2.5 such that together 
with a typical local phonon frequency of the order of 
luq = 0.1 we have a reasonable bipolaron binding energy 
Ebp = o?hu>o which pins the chemical potential at about 
half the renormalized Bosonic level = Ag - hujQa 2 . 
We choose As to lie close to the band center such that 
the number of electrons is slightly below half-filling (typi- 
cally around 0.75 per site, which is the physically relevant 
regime of concentrations). For larger binding energies 
the bipolaronic level would drop below the band of the 
electrons leading to a situation of bipolaronic supercon- 
ductivity, which is clearly not realized in HT C SC since 
they definitely show a Fermi surface. 

The idea behind applying the Boson-Fermion scenario 
to HT C SC is that we are confronted with inhomogeneous 
systems consisting of highly polarizable substructures on 
which localized bipolarons are formed. These local sub- 
structures are embedded in the rest of the lattice jj| 
which is occupied by electrons having a large either hole 
or electron-like Fermi surface depending on doping. 
In such a two-component scenario the electrons scatter in 
a resonant fashion in and out of the Bosonic bipolaronic 
states. It is that which is at the origin of the opening 
of the pseudogap in the normal state of these materials, 
driven by a precursor of electron pairing [B , rather than 
magnetic interactions Q. Generalizing this scenario in 
the way described above provides a mechanism by which 
the electrons acquire polaronic features (which, unlike 
for Bosons, are not of intrinsic nature) via the charge ex- 
change term. This term thus not only controls the open- 
ing of the pseudogap as in the original BFM but also the 
appearance of the strong incoherent contributions to the 
electron spectrum arising from phonon shake-off effects. 
Given the two-subsystem picture on which the Boson- 
Fcrmion model is based, doping leads primarily to the 
creation of localized bipolarons which beyond a certain 
critical concentration are exchanged with the itinerant 
electrons. For a system such as, for instance, YBCO the 
number tib = (bfbi) of doping induced bipolarons (ap- 
proximately given by half the number of dopant (1) 
ions in the chains) varies between and 0.5 per effec- 



tive site and the number of Fermions np = ^2 a (cf a Cicr) is 
equal to 1 if the Boson-Fermion exchange coupling were 
absent. We thus obtain a total number of charge carriers 
ntot — np + 2ns close to 2 for optimally doped systems. 
We should however remember that in real systems doping 
not only changes n to t but also the relative occupancy of 
Fermions and Bosons, which seems to be the most impor- 
tant effect in the doping mechanism of these materials, 
achieved in our model as soon as v is different from zero. 

We shall in the following solve the generalized BFM in 
the atomic limit (i.e., putting the second term in Eq.(2) 
equal to zero) for a grand canonical ensemble. In this 
case the eigenstates of the Hamiltonian are 

|0,f> = | 0)®|0)® |$(-X-)>, 

|i,j> = IT>®|o)<8|$po>, 

|2,l) = |i)(g)|0)(g)|$(X)) J 
|3,0=u;,+ |U>®|0)®|$(X)>„,, + 

+«i,+|o) ® |i) ® 

|4,i)=«i,-| U> ® |0) ® 

+«j,_|o)®|i)®|$pO) Vli _ 

|5,i) = |T)®|l)®|*(*-*o))i 

|M> = | i)®|i)® \$(x-x Q )), 

|7,0 = |TI)®|1)®|$(X-Xo))i , (3) 

where \a) denotes a site occupied by an electron with spin 
a and | f|) a site occupied by a pair of electrons with spin 
up and down. |0) and |1) denote a site unoccupied and, 
respectively, occupied by a Boson. \§{X))i denotes the 
Z-th excited oscillator state and |$(A — Xo))i = (a + — 
a) l /Vl\ exp(a(a— a + ))\$(x))o the i-th excited shifted os- 
cillator state. These two sets of oscillator states are suf- 
ficient to describe all the states listed in Eq.(3) except 
for the states |3, 1) and |4, 1) for which the corresponding 
oscillator states are given by \$(X)) Ul± and \$>(X)) Vl ± . 
The latter are determined by numerical diagonalization 
by expanding them in a set of excited harmonic oscilla- 
tor states in the form ui,±\$(X)) ul>± = T,n u L±\^( x ))n 
fwdvi,±\$(X)) vti± = J2n v "±\®( x ))n- For the regime of 
coupling parameters which we are interested in we take 
into account up to 50 phonon states, i.e., n < 50. It is 
the states |3, 1) and |4, 1) which describe the transfer of 
polaronic features from the localized bipolarons to the 
conduction electrons when Boson-Fermion exchange pro- 
cesses take place. Since photoemission only couples to 
the electrons, it is via this transfer of polaronic features 
to the intrinsically non-polaronic electrons that photoe- 
mission spectra show features which are reminiscent of 
polaronic quasi-particles. These temporarily fluctuating 
local lattice deformations described by the corresponding 
oscillator wave functions \$(X)) Ul ± and \$(X)) VI ± are 
manifest in the pair distribution function (PDF) 

9(x) = -^Y, ex v(-0E(n,l))(n,l\S(x)\n,l) . (4) 
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Here Z = YL=o Em e~^ m ^ denotes the partition 
function, with E(m, I) being the eigenvalues of the eigen- 
states listed above, given by: E(0,l) = IHloq, E(l,l) = 
E(2,l) = e + lhu 0l E{3,1) = e, i+J E(A,l) = £,,_, 
J5(5, = £7(6, = £ Q + E + lhu) - e BP and E(7, 1) = 
2eQ+Eo+lhujQ — eBP, witheo = D — fi and Eq = As — 2/j. 

In order to investigate the various physical quantities 
on the basis of this single-site generalized BFM we must 
choose As in a way to guarantee the conditions set out 
above, that is, a concentration of electrons np — 0.75 
(corresponding to a hole concentration of ~ 0.25) for a 
total concentration of particles n to t = 2. In order to 
achieve these conditions we put the bare bipolaronic level 
A b above the bare electronic energy level D such that the 
bipolaronic levelshift ebp brings this level down slightly 
below the bare electronic level. We adjust the precise po- 
sition of this level by putting As = 2D + hoj^a 2 — 5 As, 
with 8Ab = 0.025, chosen in order to give rip — 0.75. 



We have Ippsi^) — lE{u)np(uj), where np(uj) denotes 
the Fermi distribution function and Ie(w) the emission 
part of the total one-particle Fermionic spectral function 

J(a;) = i ]T ( e -/^(m,0 +e -/3£(m'/)) 

| (m', l'\ Cl \m, I) \ 2 6(uj - E(m, I) + E(m', I')) 

\2( —0(ea+m%w o ) 



:Z F 5{u-e Q ) + ± Yl Ks\ 2 (e 
+ e - /tel * a )£(u> + so + rrihojQ - ei i8 ) 

- E 



l,m,s— ± 



n< * — ' \n 

n'=0 



(I - n'\) 



-lhw + e B p) ■ (5) 



T=0.001 
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Here Z F = Ul + e - pe °){l + e -^ eo+E °- 



FIG. 1. The pair distribution function for several temper- 
atures T in units of D, the highest T corresponding to the 
most pronounced double-peak structure. The displacement is 
measured by the dimensionless parameter £ = X-\J Mwo/ft- 

Given this choice of parameters, we obtain a PDF 
(illustrated in Fig.l) showing a double-peak structure 
which merges into a single-peak structure as the tem- 
perature is lowered below the characteristic temperature 
T*, at which, as we shall see below, the pseudogap opens 
up. The two peak positions characterize the two deforma- 
tions of the local lattice environment where a given site 
is alternatively occupied by a pair of electrons or by a 
bipolaron. Recent EXAFS [§, XANES @ and pulsed 
neutron scattering |)| experiments give some indications 
for such dynamical local lattice fluctuations. 

Let us now embark on the evaluation of the inten- 
sity of the photoemission spectrum Ipes(u) from a sin- 
gle site Boson-Fermion system - tantamount to the an- 
gle integrated rather than angle resolved photoemission 
spectroscopy when neglecting the effect of the dynamical 
mean field coming from the itinerancy of the electrons. 



represents the spectral weight of the non-bonding con- 
tributions which accounts for the coherent part of the 
photoemission spectrum, unaffected by any coupling to 
the Bosons and hence to the phonons (jib(w) denotes the 
Bose distribution function). The second and third con- 
tribution to the spectral function account for the 
incoherent part of the spectrum. We illustrate in Fig. 2 
the photoemission spectral intensity IpEsi^) f° r differ- 
ent temperatures (in units of D) . For high temperatures 
(T ~ 0.06) we observe a very much broadened spectral 
function which in shape comes close to that of a typical 
Fermi liquid. Upon lowering the temperature this spec- 
tral function starts exhibiting a pseudogap and at the 
same time a broad incoherent contribution (coming from 
the second and third term of the expression for I(tv) in 
Eq.(5)) emerges. The incoherent part of the spectrum 
extends over a region in energy which is of the order of 
the half band width (~ 0.5eV) and is practically tem- 
perature independent at low temperatures, which seems 
to be confirmed experimentally p0| . The closing up of 
the pseudogap (measured as the difference in energy be- 
tween the chemical potential at lu — and the midpoint 
of the leading edge of the photoemission spectrum) as 
we increase the temperature is illustrated in the inset 
of Fig. 2. The pseudogap has a zero temperature limit 
of 0.085D ~ AOmeV and closes up at a characteristic 
temperature T* ~ 0.061? ~ 350 If, which are reasonable 
numbers. We should add that the chemical potential 
for temperatures below T* turns out to be practically 
temperature independent. In order to illustrate the clos- 
ing up of the pseudogap as the temperature approaches 
T* we plot in Fig. 3 the density of states I(oj) for dif- 
ferent temperatures. We clearly notice a strongly non- 
symmetric bias near the chemical potential (lu = 0) which 
seems to be verified in tunneling experiments H]. 
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FIG. 2. The photoemission spectral intensity Ipbs(uj), 
taking into account experimental broadening of Aw = 0.05, 
for various temperatures T, the lowest T corresponding to 
the curves with the largest pseudogap. The inset shows the 
temperature behaviour of the pseudogap. 

The work reported in this Letter relates the temper- 
ature dependence of the pseudogap to that of the inco- 
herent part of the quasi-particle spectrum. The opening 
of the pseudogap, being associated to resonant exchange 
tunneling between intrinsically unpaired electrons and 
electrons paired up in bipolaronic states, is thus driven 
by a metal-insulator cross-over rather than by supercon- 
ducting fluctuations and thus can open up well above 
the onset of the superconducting phase. The broad in- 
coherent part of this spectrum is attributed to phonon 
shake-off effects arising from the polaronic character of 
the electrons in the metallic layers, transmitted to them 
via their resonant scattering into bipolaronic states. The 
temporarily fluctuating local lattice deformations which 
are caused in this process are expected to show up in 
a characteristic double-peak feature of the pair distribu- 
tion function (measureable by EXAFS) and should test 
whether the incoherent ARPES background is of pola- 
ronic origin or not. Our approach is based on an atomic 
limit calculation of the generalized Boson- Fermion model 
which is solved exactly by numerical means. The results 
obtained are not expected to change qualitatively when 
taking into account the itinerancy of the electrons. This 
will only introduce possible asymmetries in the Brillouin 
zone coming from asymmetric coupling v in a more mi- 
croscopic model (in accordance with the hypothesis of 
strongly hybridized plane and out of plane states in cer- 
tain parts of the Brilloun zone, as suggested by LDA 
calculations [jl2|) and affect the quasi-particle structure 
close to the Fermi energy. For this energy regime our 
self-consistent studies (J] on the original BFM, fully tak- 
ing into account electron itinerancy but neglecting any 
coupling to the phonons, reproduces more faithfully the 
quasi particle structure, but with an incoherent contri- 
bution which is typically of the order of 0.1 D. This is 



one order of magnitude less than the incoherent contribu- 
tions due to phonon shake-off reported in this Letter. We 
therefore can safely treat this problem within the approx- 
imation scheme presented here. Our results are moreover 
robust in the sense that they hold for different total con- 
centrations between 1.5 and 2, as long as we enforce the 
condition that Uf is between 0.7 and 1. 
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FIG. 3. Evolution of the electron DOS I(ui) as a function 
of temperature T, the lowest T corresponding to the deepest 
pseudogap (experimental broadening Aw = 0.05). 
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